Growth and phase velocity of self-modulated beam-driven plasma waves 
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A long, relativistic charged particle beam propagating in a plasma is subject to the self-modulation 
instability. This instability is analyzed and the growth rate is calculated, including the phase 
relation. The phase velocity of the accelerating field is shown to be significantly less than the drive 
beam velocity. These results indicate that the energy gain of a plasma accelerator driven by a self- 
modulated beam will be severely limited by dephasing. In the long-beam, strongly-coupled regime, 
dephasing is reached in less than four e-foldings, independent of beam-plasma parameters. 



PACS numbers: 52.40.Mj, 52.35.-g 

Plasma-based accelerators have attracted considerable 
attention owing to the ultrahigh field gradients sustain- 
able in an electron plasma wave, enabling compact ac- 
celerators. The electric field amplitude of the electron 
plasma wave (space-charge oscillation) is on the order 
of Eq = cm e uip/e, or £7rj[V/m] ~ 96y^o|cm-^], where 
uj p = {^Kn^e 2 /me) 1 / 2 is the electron plasma frequency, 
no is the ambient electron number density, m e and e are 
the electron rest mass and charge, respectively, and c is 
the speed of light in vacuum. This field amplitude can be 
several orders of magnitude greater than conventional ac- 
celerators. Electron plasma waves with relativistic phase 
velocities may be excited by the nonlinear ponderomotive 
force of an intense laser [1] or the space-charge force of a 
charged particle beam, i.e., a plasma wakefield accelera- 
tor (PWFA) [2, 3]. In 2006, high quality 1 GeV electron 
beams were produced using using 40 TW laser pulses in 
cm-scale plasmas [4]. In 2007, a 42 GeV electron beam 
in a meter-long plasma was used to double the energy 
of a small fraction of electrons on the beam tail by the 
plasma wave excited by the beam head [5] . These exper- 
imental successes have resulted in further interest in the 
development of plasma-based acceleration as a basis for 
future linear colliders [6, 7]. 

It has recently been proposed to drive a plasma ac- 
celerator with a highly relativistic proton beam, such 
as those available at CERN (European Organization for 
Nuclear Research) [8, 9]. In general, exciting plasma 
waves requires a drive beam density profile with fre- 
quency components at the plasma frequency, i.e., a beam 
density longitudinal scale length on the order of the 
plasma wavelength X p = 2Tr/k p = 2ttc/uj p , or A p [/Km] = 
3.3 x 10 10 / ^/n[cm~ 3 ]. Compact, high-gradient acceler- 
ators require high plasma density, and therefore require 
short drive beams, e.g., X p ~ 100 for n ~ 10 17 
cm~ 3 . Generating short proton beams (or proton beams 
with spatial structure at X p ) is challenging, and it has 
been proposed to rely on a beam-plasma instability to 
modulate the beam at X p , driving a large amplitude 
plasma wave [10]. The self- modulation of the beam oc- 
curs through coupling of the transverse wakefield with 



the beam radius evolution. Periodic regions of focusing 
and defocusing modulate the beam density at X p , driv- 
ing a larger plasma density modulation that further fo- 
cuses the beam periodically. This is somewhat similar to 
the self-modulation instability that occurs for long laser 
pulses [11]. For beams long compared to X p , where self- 
modulation occurs, the instability is enabled by the drive 
beam dynamics, and therefore the wakefield properties 
will be strongly affected by the drive beam dynamics. 

An important quantity characterizing the performance 
of a plasma accelerator is the phase velocity v p of the 
plasma wave. For v p < c, a highly relativistic elec- 
tron will outrun the plasma wave and phase slip from 
the accelerating to the decelerating phase region of the 
plasma wave. This limits the electron energy gain to 
AW ~ r ) 2 (E z jE^)ra e c 1 after acceleration over a dephas- 
ing length Ld ~ 1 p X p , where E z is the electric field ampli- 
tude of the plasma wave and j p = (1 — v^/c 2 )^ 1 / 2 . For 
a plasma accelerator driven by a short (< X p ) intense 
laser pulse, v p can be relatively low (j p ~ 10 — 100) and 
dephasing can limit the energy gain [12]. For a PWFA 
driven by a short {< X p ) highly- relativistic beam, v p can 
be sufficiently high so that dephasing is not an issue. 

In this Letter we calculate the self-modulation of par- 
ticle beams in plasma, including the properties of the ex- 
cited plasma wave. In particular, we show that the phase 
velocity of the plasma wave excited by self-modulation 
is greatly reduced from the velocity of the drive beam. 
The phase velocity is determined by the growth of the 
instability and the beam-plasma dynamics. A similar ef- 
fect occurs in self-modulated laser-driven plasma waves 
[13, 14]. Analytic solutions for the growth rate and phase 
velocity in the long beam regime are derived and com- 
pared to numerical solutions of the envelope equation for 
the particle beam. Owing to the low phase velocity of 
the plasma wave, the maximum energy gain in such a 
self-modulated beam-driven accelerator will be severely 
limited by dephasing. 

The wake generated by a relativistic particle beam 
driver moving through a plasma can be calculated using 
the cold plasma fluid and Maxwell equations. Here we 
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consider a drive beam consisting of particles with charge 
=pe and mass M b . In the linear wake regime, the nor- 
malized density perturbation Sn/n = (n — n a )/n < 1 
driven by a beam with density n b < no is 



(d 2 + k 2 ) 5n/n = ^k 2 n b /n , 



(1) 



where the =F corresponds to a negatively/positively 
charged particle beam. A highly relativistic beam is as- 
sumed with Lorentz factor 7 = (1 — /J 2 ) -1 / 2 » 1, and the 
quasi-static approximation is taken such that the plasma 
fluid quantities are functions of the co-moving variable 
Q = z — fi b t. The beam-driven longitudinal electric field 
E z and transverse fields E r and Bg are [15] 

(Vi - k 2 p ) E z /E = -k p d c 8n/n , (2) 
(Vi - kl) (E r - B e )/E = -k p d r 6n/n Q . (3) 

The transverse beam-driven wakefield Eq. (3) is coupled 
to the envelope equation for the beam [16] 



d 2 R 



I TY~I 

T-5TF{k P r{E r -B e )/E Q ), (4) 



dz 2 4j 2 R 3 7 i? M b 

where R = (r 2 ) 1 / 2 is the rms beam size, e n = 
"f[(r 2 )((dr/dz) 2 ) - (rdr / 'dz) 2 ] 1 / 2 /2 is the normalized 
transverse emittance in cylindrical geometry, and the 
brackets indicate an average over the transverse beam 
distribution. 

For simplicity, in the following we consider a beam with 
a flat-top radial profile, n b = [n b or 2 /r 2 }f(()<d(r - r b ), 
where / is the normalized longitudinal profile, O is the 
Heaviside function, r b ((, z) is the beam radius, and r bn = 
r b ((,z — 0) is the initial beam radius. For a flat-top 
radial profile, Eqs. (I) and (3) have the solution 

(E r - Bg)/E = ±{n w /n )k 2 p rl h(k p r) 

x / C <sin[fc p (C-C')]/(C')^i(fcp^(C'))A 6 (C), (5) 
J 00 

for r <r bl assuming the initial radius r bn is independent 
of (. Here I m and K m are the modified Bessel functions. 
Using Eqs. (4) and (5), the envelope equation for the 
beam radius r b ((, z) — \f2R at any slice ( is 

d 2 r b e 2 4k 2 r 2 b0 I 2 (k p r b ) 



dz 2 j 2 rf 



x / d£sin[k p ((-<;')}f(<;')Ki(k p r b (C))/r b (<;'), (6) 
J 00 

where k 2 = Ann^e 2 /M b is plasma wavenumber of the 
beam. Equation (6) describes the coupled beam evolu- 
tion and wakefield excitation. 

Consider a long beam compared to the plasma wave- 
length, where the variation in the longitudinal profile 
may be neglected /(C) — 1, propagating in a plasma 



with a perturbation at the plasma wavelength. In the 
following we will decompose the plasma and beam quan- 
tities such that Q = Q + Q\, where the '0' subscripts 
indicate the long beam solution and the 'I ' subscripts in- 
dicate the perturbation. The initial perturbation or in- 
stability seed may be due to excitation of a plasma wave 
from the head of the beam, fluctuations in the beam or 
plasma, or from a plasma wave externally excited (e.g., 
by a short-pulse laser). Using Eqs. (I)-(4), the unper- 
turbed long beam solution is (5n)o = —n b , (E z )o = 0, 
(E r -B e ) = ±Eo{n bo /no)k p r K 1 (k p ro)Ii(k p r)(r b o/r ) 2 
for r < r , and the beam radius evolves as 



d 2 r 
dz 2 



7 2 r 3 



+ 



ifi(fc p r )/ 2 (fcpro) = 0. 



7Vo 



(7) 



We will assume that the beam is initially in the long beam 
equilibrium r b0 = r 
r cq is given by e 2 fc p = 



To 



< oq , such that d 2 r /dz 2 = 0, and 
47fc 2 r3 q X 1 (fc J ,r cq )/ 2 (fcp?' cq ). For a 
narrow beam, k p r eq <§C 1, the equilibrium beam radius is 
r cq = [2e 2 Jk 2 bl Y'\ 

Assuming a small perturbation about this equilibrium, 
rb = r a + n with |ri/r | <C 1, and expanding Eq. (6) 
yields the evolution of the beam radius perturbation 



df 2+4K )ri 



with the constants 
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f C dC 

J OO 



sin(C-C , )r 1 (C , ) ) (8) 



k 2 = 2K 1 (k p r ) 



.h{k p r ) ' 



(9) 



and v = 4:l2(k p ro)K2(k p r ), and the normalized vari- 
ables (■ = fcpC and z — k b z/ (2"f) 1 ' 2 . In the limit of a 
narrow beam k p rQ <C 1, v ~ 1 — (fc p ro) 2 /6, and k 2 ~ 
l + (fcp^o) 2 [C 7 -l/4 + ln(fc p r /2)]/2, where C 7 ~ 0.577 is 
the Euler-Mascheroni constant. Equation (8) may be an- 
alyzed in several regimes. The most relevant regime for 
plasma accelerators based on self-modulated drive beams 
is the strongly-coupled (or long-beam, early-time) regime 
valid for ( » z. 

Applying the linear plasma wave operator to Eq. (8) 
yields 



(d? + I)(df + 4K 2 )r 1 



2v ri . 



(10) 



Consider a slowly varying envelope, such that n = 
fexp(ik p ()/2 + c.c. with |c^f| <§C \k p f\, and assume 
the strongly- coupled regime where the growth length of 
the instability is short compared to 7 1 / 2 £; fc ~ 1 , such that 
\dzf\ ^> 2n\r\. In this regime Eq. (10) becomes 



0. 



(11) 



which describes the evolution of the slowly varying ampli- 
tude of the beam radius perturbation and may be solved 
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FIG. 1. (Color online) Beam radius modulation rt/ro vs fc p £ 
with beam-plasma parameters rib /no = 0.008, 7 = 107, and 
k p ro — 1 (and r&o = ro = r cq ), obtained from numerical 
solution of Eq. (6), at k p z = 8000 (red curve) and k v z — 
9500 (blue curve). Dashed curves are the envelope of the 
asymptotic linear solution Eq. (13). 



using standard Laplace transform techniques. With the 
initial conditions f(z,( = 0) = SrQ(z), f(z — 0, C) = Sr. 
and d z f(z = 0,C) = 0, the solution to Eq. (11) can be 
expressed as 



f/8r 



E 



2\n 



nl(2n)\ 



F 2 (;{l/2,l};z^|C|£ 2 /4), (12) 



where q F p is the generalized hypergeometric function. 
The solution to Eq. (11) may also be evaluated asymp- 
totically and has the form 



6rj^^N- 1 / 2 e N cos (k p ( + N/ V3 - n/u) 



where the number of e-foldings is 



4 V n M bl pM 



1/3 



(13) 



(14) 



Note that growth Eq. (14) [and the beam envelope equa- 
tion, Eq. (6)] differ from that found in Ref. [10]. 

Figure 1 shows the beam radius modulation r b /ro = 
1 + rj versus k p (, after propagating k p z — 8000 (red 
curve) and k p z — 9500 (blue curve), obtained from nu- 
merical solution of Eq. (6) for a beam initially in equi- 
librium r;,o = ro = r oq with beam-plasma parameters 
nb/riQ — 0.008, 7 = 107, and k p r = 1. The dashed 
curves are the envelope of the linear asymptotic solu- 
tion Eq. (13). Figure 1 shows the growth versus distance 
behind the head of the beam (at k p ( = 0) and versus 
propagation distance. Also shown is the shift in phase of 
the modulation versus propagation distance, resulting in 
a reduced phase velocity, as discussed below. 



The above solution Eq. (12) assumed \k p r\ ^> \d^f\, or 
1 3> \k~ 1 (dc i N)\. This condition may be expressed as 



1 > 



33/2 
2 6 



k p 7 



(15) 



or ( ^> z, which will be satisfied for long beams suffi- 
ciently early in the beam propagation. It was also as- 
sumed that \d&r\ 3> 2/e|r|, which is satisfied provided 
Eq. (15) is satisfied. The above analysis is also based on 
linear theory, and nonlinear effects (i.e., when r*i ~ ro or 
E z ~ Eq) may saturate the instability. 

The beam radius perturbation r*i = f exp(ifc p £)/2 + 
c.c. modulates beam density n b ~ n b (ro)(l — 2ri/ro)- 
This beam density modulation drives a modulation in 
the electron plasma density nexp(ifc p £)/2 + c.c, via 
Eq. (1), i.e., d^n ~ ^ik p n b (ro)r / r$. The plasma density 
modulation drives the accelerating wakefield E z /Eq = 
E z exp(ifc p C)/2 + c.c, via Eq. (2), i.e., (Vi - k%)d ( E z = 
z Fkp[nt > (ro)/no]f/rQ. For the same initial conditions as 
above, the series solution for the accelerating wakefield 
in the long-beam regime is 



E, 



TH R {r,r ) \Q V 7 — , lWo 



(16) 



and the sum may be expressed as the hypergeometric 
function oF 2 (;{l/2,2};w|C|z 2 /4). Here H R (r,r ) = 1 - 
k p r a Ki(k p ro)Io(k P r ) for r < r o and k p r I 1 (k p r Q )K (k p r) 
for r > ro. In the asymptotic limit, E z /E z (z — 0) ~ 
3 7/4 (32n)- 1 / 2 N~ 3 / 2 exp(TV) cos^), where the number of 
e-foldings of growth of the accelerating wake is given by 
Eq. (14) and the phase is 



V> = k pC - T + 
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1/3 



(17) 



The phase velocity of the accelerating wake is given by 
f3 p = -d t ^/d z ip = d c ijj/(d c + d z )^ ~ 1 - drf/drf. In 
this regime, i.e., satisfying Eq. (15), f3 p ~ 1 — k~ x d z ii. 
Using the phase Eq. (17), the phase velocity is (3 p = 
1 - {2/3 3 / 2 )(N/k p z). The phase velocity of the self- 
modulated beam-driven wakefield is less than the beam 
velocity /3 b ~ 1, varies along the beam £ and during 
propagation z. Asymptotically, the Lorentz factor of the 
phase velocity is 



7 P 



/ jn M b z 
\vn bo m e |(| 



1/6 



(18) 



in the strongly-coupled, long-beam regime. Note that, 
behind the modulated beam the phase velocity is given 
by Eq. (18) with \C\ = L b , where L b is the bunch length. 

Figure 2 shows the normalized Lorentz factor 
of the phase velocity of the accelerating wakefield 
7p[i/(fcf,/fcp) 2 |C|/7] -1 / 4 versus normalized propagation 
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FIG. 2. (Color online) Normalized phase velocity of accelerat- 
ing wakefield Jp[v(kb/k p ) 2 \£\/'y]~ 1 / 4 vs. normalized propaga- 
tion distance (^|£|) 1,/2 z in the long-beam regime: solid (black) 
curve is the series solution Eq. (16), dashed (red) curve is the 
asymptotic solution Eq. (18), and (blue) dots are from the 
numerical solution of the envelope equation Eq. (6). 



distance (^l^l) 1 / 2 ^. The solid curve in Fig. 2 is 
obtained from the series solution Eq. (16), (3 p = 
1 - k- l d z [a^ct&n{^E z /^.E z )], the dashed curve is the 
asymptotic solution Eq. (18), and the dots are from 
the numerical solution (with the parameters 7 = 107, 
n-bo/no = 0.008, k p r^ = 1, and k p L b — 600) of the en- 
velope equation Eq. (6). Figure 2 indicates that there is 
a minimum phase velocity. The minimum phase velocity 
can be estimated by using the series solution Eq. (16). 
The minimum phase velocity occurs at (^l^l) 1 / 2 ^ ~ 1.72, 
with 



1.06 



/ 7n M b 
\vn w m e k p \C,\ 



1/4 



(19) 



As shown in Fig. 2, after reaching 7 m i n , the phase ve- 
locity grows slowly as the beam propagates 7 P oc z 1 / 6 
[cf. Eq. (18)]. For example, consider a wake driven 
by a 100 GeV proton beam (7 = 107), with tq = 
180 /mi, L b = 10 cm, and 10 11 particles. Operating 
at no — 10 15 cm~ 3 , corresponds to n b o/no = 0.008, 
E = 3 GV/m, k p r = 1.0, k p L b = 600, and v = 0.88. 
For this example, 7 m i n — 15 behind the drive beam after 
z ~ 8.5 cm (i.e., w 500 A p ) of propagation. 

With the phase velocity of the self-modulated wake 
determined, the dephasing length may be calculated. 
For a linear wake, the dephasing length is the propa- 
gation distance required for an ultra-relativistic particle 
/3b ~ 1 to slip Ap/4 (or a wake phase of 7r/2) with re- 
spect to the plasma wave. Assuming the phase velocity 
is well-approximated by the asymptotic solution in the 
strongly-coupled regime Eq. (18), the dephasing length 
is L d = (27r/3) 3 / 2 (^fc 2 fc p |C l |/7)" 1/2 . Including the early 



time response via Eq. (16), the dephasing length is 

L d ~ 4.9 (i/fcgfcplG I/7) ~ 1/2 , (20) 

where Q is the injection position of the witness bunch 
(e.g., initially at a peak of the accelerating field). For a 
witness bunch injected behind the drive beam, Q = L b . 
This reduced dephasing length will greatly limit the en- 
ergy gain of a witness electron beam trailing the drive 
bunch. For example the number of e-foldings of the self- 
modulational instability that have occurred at the de- 
phasing length Eq. (20) is N(z = L d , &) ~ 3.8. Note that 
the number of e-foldings at a dephasing length N(z = Lg) 
is independent of injection location and the beam-plasma 
parameters. Improved efficiency may be possible by ta- 
pering the plasma density, i.e., increasing the background 
plasma density to reduce the plasma wavelength, thereby 
increasing the phase velocity [17], although variation of 
the density may affect the instability growth. Alterna- 
tively the accelerator may use a staged approach, where 
a long plasma region self-modulates the drive beam, fol- 
lowed by a second stage where a witness bunch would be 
injected following the modulated drive beam. Such a two- 
staged approach could potentially be limited by the hose 
(or transverse two-stream) instability [18], which grows 
in the long beam limit with a comparable growth rate 
~ N. This implies that to drive large amplitude acceler- 
ating fields via the self-modulational instability without 
hosing requires strongly seeding the instability. One pos- 
sibility is to use a beam with a fast rise in the current 
profile [10]. Another possibility to seed the modulation 
is via an intense short-pulse laser. 

The long-beam, early-time regime described above will 
be valid for z <C Q. After sufficiently long propagation 
distances, or for sufficiently short beams, the instability 
may enter a weakly-coupled regime where the instability 
growth length is long compared to 7 1 / 2 jk b . The insta- 
bility will transition to the weakly-coupled regime after 
a propagation distance approximately z ~ Q, or, using 
Eq. (14), after approximately N ~ k p (. For long beams 
k p Q 3> 1, nonlinear effects will typically appear before 
the instability enters the weakly-coupled regime. 

In this Letter we have calculated the beam self- 
modulation instability growth rate, in the long-beam 
regime, including the phase dependence. The phase ve- 
locity of the accelerating wakefield was calculated and 
shown to be significantly less than the drive beam ve- 
locity. The dephasing length was calculated, and, in the 
strongly-coupled regime, a witness beam will reach de- 
phasing in less than four e-foldings, independent of beam- 
plasma parameters. This indicates that the energy gain 
in a plasma accelerator driven by a self-modulated PWFA 
will be limited by dephasing. 
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